Spin-orbit Coupled Bose-Einstein Condensates in Spin-dependent Optical Lattices 
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We investigate the ground-state properties of spin-orbit coupled Bose-Einstein condensates in spin- 
dependent optical lattices. The competition between the spin-orbit coupling strength and the depth 
of the optical lattice leads to a rich phase diagram. Without spin-orbit coupling, the spin-dependent 
optical lattices separate the condensates into alternating spin domains with opposite magnetization 
directions. With relatively weak spin-orbit coupling, the spin domain wall is dramatically changed 
from Neel wall to Bloch wall. For sufficiently strong spin-orbit coupling, vortex chains and antivortex 
chains are excited in the spin-up and spin-down domains respectively, corresponding to the formation 
of a lattice composed of meron-pairs and antimeron-pairs in the pseudospin representation. We also 
discuss how to observe these phenomena in real experiments. 
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Introduction. — In recent years, the experimental con- 
trol on ultracold atomic gases has reached truly unprece- 
dented levels. By employing two lasers with different 
frequencies and polarizations plus a non-uniform vertical 
magnetic field, experimentalists have produced spin-orbit 
(SO) couphng, which couples the internal states and the 
orbit motion of the atoms [1-4]. The SO coupled ultra- 
cold atomic gases have attracted great interests of re- 
searchers [5-15]. It has been indicated that the interplay 
among SO coupling, interatomic interaction and exter- 
nal potential leads to rich ground-state phases, such as 
plane wave, density stripe, fractional vortex and various 
vortex lattices [16-26]. The SO coupled ultracold atomic 
gases open a new window for quantum simulation, and 
provide opportunities to study SO coupling phenomena 
in a highly controllable impurity- free environment. 

All the existing studies on SO coupled Bose-Einstein 
condensates (BECs) only refer to the case that different 
internal states of the atoms are trapped in an identical 
external potential. However, by using two counterpropa- 
gating lasers with the same frequency but different polar- 
izations, the experimentalists have been able to produce 
spin-dependent optical lattices that allow different inter- 
nal states of the atoms experience drastically different 
external potentials [27-31]. The spin-dependent optical 
lattices bring more complicated geometry of condensates 
and have potential applications in quantum computation 
[27], cooling and thermometry [28], and quantum simu- 
lation [29-32]. It is natural to ask what new structures 
can be formed due to the competition between the SO 
coupling and the spin-dependent optical lattices. 

In this Letter, we investigate the ground-state phase 
diagram of SO coupled BECs in spin-dependent optical 
lattices. In the absence of SO coupling, the ground state 
of the system is characterized by the formation of al- 
ternating spin domains. However, such a structure can 
be dramatically changed due to the SO coupling effects. 
Relatively weak SO coupling basically changes the ori- 



entation of the spins in the domain walls, causing the 
transformation from Neel wall to Bloch wall. Sufficiently 
strong SO coupling excites meron-pairs and antimeron- 
pairs in the spin-up and spin-down domains respectively 
and generates a meron-pair lattice. This is essentially 
different from the mechanism of generating a meron-pair 
lattice by bulk rotation [33]. Our findings provide a new 
way to create and manipulate topological excitations in 
SO coupled systems. 

Energy hand structure. — We consider SO coupled 
BECs confined in the combined potential of a quasi-2D 
harmonic trap and ID spin-dependent optical lattices. 
The Hamiltonian of this system is given by 
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where ^ = [^^(r), ^^(r)]^ denotes the two-component 
wave functions and is normalized as J dr^^^ = N with 
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FIG. 1: (Color online) (a) Single-particle phase diagram 
spanned by the SO coupling strength k, and the depth Vb of 
the optical lattice, (b) Many-body phase diagram spanned by 
K and Vb with the effective interaction parameter g = 6000. 
The dashed line in (b) indicates the phase boundary between 
phase I and phase II in (a) for comparison. 



FIG. 2: (Color online) Band structures induced by the competition between the SO coupling strength n and the depth Vo of 
the optical lattice corresponding to phase I (a), phase IIA (b), and phase IIB (c) of the single-particle phase diagram in Fig. 
1(a). The superposition of the momenta at the minima of the bands produces a density stripe in (a), a lattice in (b) and (c). 



the total particle number, = ^*^| and = ^*^| 
describe the particle number density of each component. 
gij = 47rh^aij/m {i^j = 1,2) represent the interatomic 
interaction strengths characterized by the s-wave scat- 
tering lengths aij and the atomic mass m. 

We consider a Rashba SO coupling Vso = —iTin{cFxdx^ 
(Jydy)^ where (Jx^y are the Pauli matrices and de- 
notes the SO coupling strength. The combined exter- 
nal potential Vi{T) = FH(r) + ^OLi(^), where Vr = 
^muj'^[{x'^ + y'^) + A^z^] is the harmonic trapping po- 
tential with A = 0Jz/^± ^ and Vqli = Vbsin^(i/x) 
and VbL2 = Vbcos^(i/x) describe the ID spin-dependent 
optical lattice potentials, which are experienced by the 
two components, respectively. Approximating the z de- 
pendence of the wave functions by the single-particle 
ground state in a harmonic potential, one can obtain 
the 2D dimensionless effective interaction parameters 
gij = 2\/2'n\Naij / ah-, with an = {muj±) [34]. 

The single-particle energy bands are critically impor- 
tant to understand the ground-state properties of the 
condensates. Without considering the harmonic trap, 
the 2D single-particle wave functions in the k-space 
obey the secular equations -^[{kx -\- 2siy) + ky]Us — 
^ {Us+i - 2Us + Us-i) + (kx + 2sv - iky) Ws = eUs 
and ^[{kx + 2svf^kl]Ws^^ (Ws+i + 2Ws + Ws-i)^ 
hhi {kx -\- 2su -\- iky) Us = sWs^ where Us = 
U {kx + 2sh'^ ky) and Ws = W {kx + 2sh'^ ky) with 
s = 0, ±1, ±2, ±00 represent the wave functions at 
the point {kx -\-2su,ky). Typically, we choose u = a^^ 
for our present discussion. By numerical exact diagonal- 
ization, we can solve the secular equations and obtain 
the energy band structure. 

We find that there exist three different kinds of energy 
band structures depending on the competition between 
the SO coupling strength and the depth Vq of the op- 
tical lattices. Fig. 1(a) presents the single-particle phase 
diagram spanned by and Vq. In phase I, the minima of 
the energy bands locate in a set of k points with ky = 
and kx e Ki = ±3^/, ±5i/, ...} [See Fig. 2(a)]. In 
phase II, the minima of the energy bands locate in a set 
of k points with ky = ±S {O < S < ^) , and kx e Ki for 
phase IIA {kx ^ K2 = {0, ±2i/, ±4z/, ±6i^, ...} for phase 



IIB) [See Figs. 2(b) and 2(c)]. 

The single-particle ground state in phase I is nonde- 
generate and can be expressed as a linear superposition 
of the plane waves with wave vectors {kx G Ki^ky = 0). 
This yields alternating spin domains with opposite mag- 
netization directions (density stripe). In both phase IIA 
and phase IIB, the single-particle ground state is double- 
degenerate. Each degenerate state can be expressed as a 
linear superposition of the plane waves with wave vectors 
{kx ^ K^ky = S) or {kx ^ K^ky = —6), where K = Ki 
for phase IIA and K = K2 for phase IIB. For an arbi- 
trary nonzero superposition of the two degenerate states, 
lattice will be formed as the single-particle ground state. 

Phase diagram. — By using the imaginary time evolu- 
tion method, we can solve Eq. (1) to obtain the many- 
body ground state. Considering that the spin-exchange 
integrations are very weak in typical experiments, we just 
discuss the case that gij = g. The many-body phase di- 
agram spanned by and Vq with g = 6000 is presented 
in Fig. 1(b). We find that the competition between the 
SO coupling strength and the depth of the optical lattice 
leads to three distinct phases — density stripe, triangular 
vortex lattice and rectangular vortex lattice. 

In the density stripe phase, the spin-up and spin-down 
components are arranged alternately and form alternat- 
ing spin domains [See Figs. 3(al-a6)]. Comparing the 
phase diagrams in Figs. 1(a) and 1(b), we find that the 
interaction has no significant influence on the phase re- 
gion of the density stripe. 

In the vortex lattice phases, both the triangular 
and rectangular lattices are composed of alternately ar- 
ranged vortex and antivortex chains, which are excited 
in the spin-up and spin-down domains respectively [See 
Figs. 3(bl-b6) and 3(cl-c6)]. The only difference is that 
the vortices of the neighboring chains are staggered for 
the triangular lattice, but are parallel for the rectan- 
gular lattice [See Figs. 3(b3) and 3(c3)]. These two 
different arrangements of vortices correspond to odd- 
parity and even-parity distributions of the particles in 
kx direction of the k-space. [See Figs. 3(b6) and 3(c6)]. 
These correspond to the single particle band structures 
described in Figs. 2(c) and (b), where the minima of 
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FIG. 3: (Color online) Ground state as a function of the SO coupling strength n and the depth Vb of the optical lattice with 
K = 4ahOJ±, Vb = 40?i(jj± (al-a6), k = AauOJ^^ Vb = 2()?iuj± (bl-b6), n = 4ahOJ±, Vb = 15?i(jJ± (c1-c6). The effective interaction 
parameter is fixed at ^ = 6000. The spin- up, spin-down, and total density profiles are shown in (al, bl, cl), (a2, b2, c2), and 
(a3, b3, c3), respectively. The phases of the spin-up and spin-down wave functions, with values ranging from —27v to 27v (blue 
to red), are shown in (a4, b4, c4) and (a5, b5, c5). The momentum distributions are depict in (a6, b6, c6). 



the bands also show odd-parity and even-parity distribu- 
tions respectively, although their phase regions are not 
consistent due to the influence of the interatomic inter- 
actions [See Figs. 1(a) and 1(b)]. As discussed above, 
the single-particle ground state in phase II is double- 
degenerate. From Figs. 3(b6) and 3(c6), we can see that 
the interaction removes the degeneracy and chooses an 
equal weighted linear superposition of the two degener- 
ate states as the many-body ground state. 

The alternating arrangement of the vortex and an- 
tivortex chains leads to alternating-direction plane waves, 
which propagating on two sides of each chain [See 
Figs. 3(b4,b5) and 3(c4,c5)]. The vortex line density 
riv and the wave number of the plane waves ky satisfy 
fi = Numerical simulations indicate that for a e^iven 
SO couphng strength as the lattice depth Vq increases 
from 0, ky gradually decreases from ^ and eventually 
becomes on the boundary of the vortex lattice phase. 
This implies that by adjusting the depth of the optical 
lattice, one can continuously control the vortex line den- 
sity from to 

Spin domain wall. — The separation between the spin- 
up and spin-down domains is not sharp, but requires the 
spin density vector varying gradually across the oppo- 
site domains and forming a spin domain wall [35] . There 
are two basic types of domain walls, Neel wall and Bloch 



wall. In the Neel wall spin flip occurs in a plane, while 
in the Bloch wall the spin flip occurs by tracing a he- 
lix [See Figs. 4(bl) and 4(b2)]. An intriguing finding 
of the present work is that the SO coupling dramati- 
cally changed the domain wall from Neel wall to Bloch 
wall. Figs. 4(al) and 4(a2) show the spin density vector 
S = without and with SO coupling. We can see 

that in the absence of SO coupling, the spin density vec- 
tor across the opposite domains forms a Neel wall, while 
in the presence of SO coupling it forms a Bloch wall. 

This phenomenon can be understood as follows. The 
direction of the spin flip in the domain wall only depends 
on the relative phase, and can be represented by an az- 
imuthal angle a = diTctdin{Sy/ Sx) = 0^ — 9^^ where 0^ and 
0^ are the phases of the wave functions. When the SO 
coupling is absent, there is a constant phase difference 
or TT [See the solid line of Fig. 4(c)], so the spin in the 
walls just flips along the x-direction and forms Neel walls. 
When the SO coupling is present, the phase difference is 
changed into [See the dashed line of Fig. 4(c)], so 
the spin in the walls just flips along the ^-direction and 
forms Bloch walls. 

Meron-pair lattice. — The regular triangular or rectan- 
gular vortex lattice obtained in Fig. 3 can be equivalently 
described by the spin density vector S in the pseudospin 
representation. Fig. 5(a) presents the vectorial plots of S 
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FIG. 4: (Color online) (a) The vectorial plots of the pseu- 
dospin S projected onto the x-y plane with g = 6000, Vb = 
40huj±, and k = (al), k = Aauoo^ (a2). The colors ranging 
from blue to red describe the values of the axial spin Sz from 
— 1 to 1. (b) 3D renderings of Neel wall and Bloch wall, (c) 
Section views of the relative phase 0^ — 0^ along the x axis 
with HI — ^ (solid line) and hi — Aanoo^ (dashed line). 



under a pseudo-spin rotation, which corresponding to the 
state represented in Figs. 3(cl-c6), and the correspond- 
ing topological charge density q (r) = ^e*-^S • 9^8 x djS 
is plotted in Fig. 5(b). One can see that the spin tex- 
ture in Fig. 5(a) represents a lattice composed of meron 
pairs and antimeron pairs [36, 37]. Either a meron pair 
or an antimeron pair has a "circular-hyperbolic" struc- 
ture [See Figs. 5(cl) and 5(c2)], and the only difference 
is that they have exactly opposite spin orientations. The 
spatial integral of q{v) indicates that a meron pair just 
carries topological charge 1, while an antimeron pair car- 
ries topological charge —1. Previous studies indicated 
that stable meron-pair lattice can be obtained in a rotat- 
ing system [33] . Our results show that meron-pair lattice 
can also be stabilized in alternating spin domains by SO 
coupHng without rotation. 

Experimental proposal. — In real experiments, we 
can choose a two-level ^^Rb BEG system with 
|F = l,m/ = l) and |F = l,m/ = -l). About TV = 
1.7 X 10^ atoms are confined in a harmonic trap with the 
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FIG. 5: (Golor online) (a) The vectorial plots of the pseu- 
dospin S projected onto the x-y plane under a pseudo-spin 
rotation Gx —(Jz and Gz cfx- (b) The topological charge 
density g'(r). (c) The amplification of the two kinds of ele- 
ments in fa). 



trapping frequencies uo^ ^ 27r x 40 Hz and cj^ ^ 27r x 200 
Hz. It is convenient to produce spin-dependent optical 
lattices with a large lattice spacing by using a CO2 laser 
operated at a wavelength of 10.6 /im. Under typical 
experimental conditions, the s-wave scattering lengths 
aij ^ lOOttB {glb is the Bohr radius). Based on these 
experimental parameters, we can calculate that the effec- 
tive interaction parameter g ^ 6000 and the wave vector 
u = a^^, which are consistent with our present calcula- 
tion. 

For a given SO coupling strength Hi = A^Tiuj^/m., 
adjusting the lattice depth from Vb = /c^ x 80 nK to 
Vq = ks >^ 40 nK, then to Vq = ks >^ 30 nK {ks is 
the Boltzmann's constant), one can directly observe the 
phase transitions from density stripe to triangular vortex 
lattice, then to rectangle vortex lattice by monitoring in 
situ the density profile. And for a given lattice depth 
Vb = X 80 nK, adjusting the SO coupling strength 
from = to /t: = 4:^Tiuj^/m^ one may indirectly ob- 
serve the transition of the spin domain wall from Neel 
wall to Bloch wall by dual state imaging technique, which 
can spatially resolve the relative phase [38]. 

Conclusion. — In summary, we have investigated the 
ground-state phase diagram of spin-orbit coupled BECs 
in spin-dependent optical lattices. We observe novel 
spin-orbit coupling effects on alternating spin domains 
produced by the spin-dependent optical lattices. Actu- 
ally, alternating spin domains exist in many physical sys- 
tems, such as magnetic materials and condensed matter 
systems [39, 40], thus similar spin-orbit coupling effects 



5 



would be discovered in those systems. We hope that our 
findings will deepen the understanding of spin-orbit cou- 
pling phenomena and provide thoughts on engineering 
new quantum states in spin-orbit coupled systems. 
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